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ABSTRACT 

Context. Determining the metal content of low-mass members of young associations provides a tool that addresses different issues, 
such as triggered star formation or the link between the metal-rich nature of planet-host stars and the early phases of planet formation. 
The Orion complex is a well known example of possible triggered star formation and is known to host a rich variety of proto-planetary 
disks around its low-mass stars. Available metallicity measurements yield discrepant results. 

Aims. We analyzed FLAMES/UVES and Giraffe spectra of low-mass members of three groups/clusters belonging to the Orion 
association. Our goal is the homogeneous determination of the metallicity of the sample stars, which allows us to look for [Fe/H] 
differences between the three regions and for the possible presence of metal-rich stars. 

Methods. Nine members of the ONC and one star each in the A Ori cluster and OB lb subgroup were analyzed. After the veiling 
determination, we retrieved the metallicity by means of equivalent widths and/or spectral synthesis using MOOG. 
Results. We obtain an average metallicity for the ONC [Fe/H] = -0.01 ± 0.04. No metal-rich stars were detected and the dispersion 
within our sample is consistent with measurement uncertainties. The metallicity of the A Ori member is also solar, while the OB lb 
star has an [Fe/H] significantly below the ONC average. If confirmed by additional [Fe/H] determinations in the OB lb subgroup, this 
result would support the triggered star formation and the self-enrichment scenario for the Orion complex. 

Key words. Stars: pre-main sequence ~ Stars: abundances - Open Clusters and Associations: Individual: Orion - Stars: planetary 
systems 



1. Introduction 

Tfie Orion OB 1 association, located at rouglily 400 pc, is proba- 
bly tfie best known region of active star formation. It comprises 
four major subgroups of different ages, sizes, and positions (la, 
lb, Ic, Id where the Orion nebula cluster -ONC- is), together with 
' several grouping of stars and clusters, covering a wide range 
of ages and environmental conditions: extremely young groups 
of stars are present still embedded in their parental clouds, 3-5 
Myr old clusters like cr and A Ori, as well as older populations 
like the Ori la sub-association (age ~ 10 Myr). Orion is since 
long considered as a possible example of triggered star forma- 
tion (Blaauw 1991, ASI Ser 342, p. 125), supported by the in- 
creasing ages between the Id/Ic, lb, and la subgroups and by the 
different content of gas and dust of the different regions. In the 
triggered scenario, winds and SN-driven shock waves originat- 
ing from a first generation of massive stars induce the formation 
of a new group/cluster of stars. This second generation popula- 
tion is likely contaminated by the enriched ejecta of the massive 
stars of the first generation (e.g.. Reeves IT978I I and thus might 
show a different abundance pattern. The secure and homoge- 
neous determination of the chemical composition across Orion, 
and the comparison of the abundance patterns of the different 
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subgroups, hence represent a critical and independent test of this 
scenario. 

Besides providing insights into triggered star formation 
models, the determination of the metallicity of low-mass stars 
in young associations is important in the context of extra so- 
lar planets. It is now well-established that planet-host stars are 
on average more metal-rich than stars without planetary-mass 
companions (Gonzalez 119971 Gonzalez et al. 120011 Santos et 
al. '2001'; Fisher & Valenti 120051) . The mean metallicity of stars 
hosting a gas giant planet is [Fe/H]= -1-0.15 + 0.23, compared to 
the value [Fe/H]= -0.1 +0.18 for the solar neighborhood. Also, 
the frequency of giant planets around metal-rich stars ([Fe/H]> 
0.3) is fp > 20%, much higher than that of planets around so- 
lar (or under-solar) metallicity stars (fp ~ 3%). More important, 
complementary studies suggest that the high metal content is 
primordial and not due to pollution by cannibalized planetary 
bodies (Ecuvillon et al. ,2006. Gilli et al. 120061 Pinsonneault et 
al. 1200 II see however Pasquini et al. l2007l for a different view). 

Surveys of the properties of planet-host stars provide clues 
to the mechanism of planet formation based on its final product 
(i.e., the already formed planet), since the majority of field stars 
where a giant planet has been detected are older than ~ 1 Gyr 
On the other hand, circumstellar disks around pre-main sequence 
(PMS) stars, which are the birth sites of planets, are ubiquitously 
found in star-forming regions (SFRs). The question then is what 
the metallicity of low-mass stars that are likely to be forming 
planets now in SFRs and young associations might be. To defini- 
tively answer it, the metallicity needs to be derived, as much as 
possible, for a complete SFR sample. The Orion complex al- 
ready represents an excellent target, given the high frequency 
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and wide variety of proto-planetary disks around its low-mass 
members. 

In the last 50 years, hundreds of studies have investigated 
key aspects related to the Orion stellar population, as well as 
to the associated HII regions and the surrounding interstellar 
medium. However, considerably less attention has been given 
to determining its chemical composition. Available information 
on the abundance pattern in Orion and, in particular, of the ONC, 
mostly comes from studies of the ionized gas in the Orion neb- 
ula or from abundance determination in stars of early spectral- 
type. The latter, however, is affected by several uncertainties, 
for example, NLTE and rotational broadening of the spectral 
lines. Few measurements of the metallicity among lower-mass 
stars have so far been carried out, as we discus in Sect. |2] be- 
low. With this background in mind, we exploited the available 
high-resolution spectra acquired in the framework of two differ- 
ent projects (Palla et al. 12005112007 1 Sacco et al. 120081) to mea- 
sure the metallicity in a sample of ONC members and in one star 
belonging to the A On cluster. One member of the Orion OB lb 
subgroup, whose spectrum was retrieved from the ESO archive, 
was also analyzed. 

The paper is organized as follows. In Sect. 2 we give an 
overview of previous determination of metallicity in Orion mem- 
bers; in Sect. 3 we describe the sample, observations, and data 
reduction. The analysis is presented in Sect. 4 and the results and 
discussion in Sect. 5, followed by the conclusions (Sect. 6). 

2. Previous metallicity determinations in Orion 

Based on the determination of oxygen in the Orion nebula by 
Osterbrock et al. C 1992 J , it has been believed until recently that 
the metallicity of Orion is lower than that of the Sun, represent- 
ing a problem for Galactic evolution models in the solar vicin- 
ity. While the reappraisal of oxygen in the Orion nebula by the 
study of Esteban (2004) no longer supported this view, stellar 
abundances obtained by Cunha & Lambert (.1994) and Cunha 
et al. d 1 9981 ) also suggested a lower-than-solar metallicity for 
the whole Orion complex. Specifically, Cunha and collaborators 
measured Fe, O, C, N, and Si abundances in a sample of B- and 
F/G- type Orion members, primarily belonging to the subgroups 
la and Ic. They found that the metallicity ([Fe/H]) ranges be- 
tween -0.31 (Ic, G-type star) and -1-0.14 (la, B-type star); in spite 
of this rather wide interval Cunha and collaborators concluded 
that all the four subgroups la, lb, Ic, and Id (where the ONC 
is located) are characterized by the same, somewhat below solar, 
metallicity ([Fe/H] = -0.08+0.13, where the error is the Icrdevi- 
ation from the mean). F and G-type stars only would indicate an 
even lower metal content ([Fe/H]= -0.13 ± 0.13). Interestingly, 
the most metal-rich among the F/G-type stars is the only ONC 
members included in the sample. For this star, PI455, they de- 
rived a metallicity [Fe/H]= -1-0.08 + 0.15. At variance with their 
claim of a homogeneous iron content, Cunha and collaborators 
( 119921 119981 ) found evidence of star-to-star variations in O and 
Si abundances, with a few O- and Si-enhanced stars being more 
centrally located (in the Trapezium region) than the O- and Si- 
poor ones, distributed throughout the association. The authors 
suggest that the observed abundance pattern could be evidence 
of the expected self-enrichment. More recently, Simdn-Diaz et 
al. (|2006i derived oxygen abundances for a few B-type stars 
in the Trapezium finding lower values than those of Cunha & 
Lambert and, in contrast to their results, no evidence of any el- 
emental enhancement. 

Focusing on low-mass stars, seven objects in Orion Ic and 
Id were included in the study of Padgett (1 19961 1, who derived 



[Fe/H] values in the interval -0.1 to -1-0.23 dex and an av- 
erage [Fe/H] slightly above solar Of the four ONC stars in- 
cluded in the sample, two showed slightly over-solar abundances 
([Fe/H]= 0. 14 ± 0. 18 and 0.08 + 0. 17 respectively), and the other 
two stars a solar metallicity ([Fe/H]= -0.01+0.21,-0.01+0.17). 
Finally, based on three stars, Santos et al. (2008) obtained an av- 
erage for the ONC [Fe/H]= -0.13 + 0.06, again below solar. 
In summary, existing determinations of the metal content in the 
whole Orion complex reach controversial conclusions both on 
the overall metallicity -in particular on whether it is below solar 
or not- on the presence of a star-to-star scatter within a given 
subgroup (the ONC in particular), and on differences in [Fe/H] 
between the different subgroups. 



3. Sample stars and observations 

3.1. ONC 

The spectra were acquired in two observing sessions with 
FLAMES on VLT/UT2 (Pasquini et al. 2000) as part of a project 
focused on measurements of lithium abundances in very low- 
mass members of the ONC (Palla et al. (20051 l2(j07l l. In both 
runs Giraffe fibers were mostly allocated to M-type stars; in- 
stead, UVES fibers were used to observe K- and M-type stars 
with the goal of deriving their accretion rates by means of the Ca 
II infrared triplet (Flaccomio et al., in prep). The present sample 
includes all the K-type stars observed with UVES (seven stars), 
while M dwarfs were discarded since determination of metallic- 
ity in these cool stars is not feasible with standard methods/codes 
and Kurucz model atmospheres (see Sect. 14. 3.T] below). To the 
UVES sample, we added the only two K-type stars observed 
with Giraffe and for which Palla et al. (I2005I I derived a negli- 
gible veiling. The total sample hence consists of nine stars. 

The observations were carried out on 15 February (Visitor 
Mode) and 15 November (Service Mode), 2004. The UVES 
spectrograph (Dekker et al. 2000) in fiber mode provides a res- 
olution R=40000. We used the CD4 cross-disperser (needed to 
target the Ca II triplet), allowing us to cover a spectral range from 
-6700 A to 10000 A. Giraffe was used with the HR15 setup, 
yielding a resolution R~ 19000 and spectral coverage from 660.7 
to 696.5 nm. 

The sample stars are listed in Table [1] where we give the iden- 
tifier from Hillenbrand ( 119971 ), an observing log (the dates of 
the observations, total exposure time, S/N ratios), along with 
the stellar properties (I magnitude, V-I color, spectral-type, ro- 
tational period, effective temperature -Tej-, and K-band ex- 
cess), are also listed in the table. Information on photometry, 
spectral-types, effective temperatures, and excess were taken 
from Hillenbrand ( fT997] l and Hillenbrand et al ( [19981 ), while 
rotational periods come from Herbst et al. ( 120021 1 . 

Reduction of UVES data was carried out using the 
FLAMES-UVES pipeline, following the standard procedure: 
bias subtraction, flat-field correction, order extraction and wave- 
length calibration. Sky subtraction was performed outside the 
pipeline context, using the fiber allocated to the sky. Spectra 
of the same star obtained in different exposures were then co- 
added, after checking that no radial velocity variations were 
present. Giraffe data were reduced using the BLDRS software 
(Blecha & Simond [20041 -see details in PaUa et al. 2007). Final 
S/N ratios are in the range 70-125 for the UVES spectra and 
150-200 for the Giraffe ones. All UVES spectra are shown in 
Fig. [I] 
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Fig. 1. UVES spectra of the ONC sample stars in the wavelength range 6740-6780 A. Different features are marked. Note the 
shallower lines in the spectrum of star h278, for which we derived a high veiling value. 



Table 1. ONC target stars, observing log, and target properties . 



Star 


Obs. Dates 


Instrument 


Exp. time 


S/N 


I 


V-I 


Sp. Type 




A(I-K) 


P 


Vmd. 


(H97) 






(sec) 










(K) 


(mag) 


(days) 


(km/s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


223a 


15-02-04 


UVES, CD4 


3x3600 


70 


13.64 


2.42 


K5 


4395 


0.48 




25.43 ± 1.58 


268 


15-02-04 


UVES, CD4 


3x3600 


95 


12.88 


1.61 


K5-K6 


4197 


0.16 


9.81 


27.80 ± 1.02 


278 


15-02-04 


UVES, CD4 


3x3600 


70 


13.66 


1.94 


K2-K7e 


4395 


1.69 


6.76 


24.58 ± 1.49 


487 


15-02-04 


UVES, CD4 


3x3600 


105 


12.79 


1.78 


K6 


4197 


0.28 




29.27 ± 0.87 


673 


15-02-04 


UVES, CD4 


3x3600 


110 


13.07 


1.83 


K5 


4395 


0.81 


1.44 


24.53 ± 1.29 


683 


15-02-04 


UVES, CD4 


3x3600 


130 


11.76 


1.85 


K6 


4197 


0.28 


11.50 


25.73 ± 1.51 


907 


15-11-04 


UVES, CD4 


4x2600 


125 


12.77 


1.60 


K1-K4 


4775 






20.82 ± 1.24 


664 


15-02-04 


GIRAFFE 


3x3600 


150 


14.05 


2.88 


K5.5-K7 


4197 


-0.18 


7.20 


25.36 ± 1.18 


1020 


15-02-04 


GIRAFFE 


3x3600 


200 


12.65 


1.40 


K6 


4197 


-0.01 




27.60 ± 1.82 



' Identifier (1), I, V-I, spectral-types (6, 7, 8), and Teff values (9 -from spectral-types using the scale of Kenyon & Hartmann 1995) come from 
Hillenbrand J1997b : K-band excess (10) and rotational periods (11) were taken from Hillenbrand et al. t l998 j and Herbst et al. C002J . respectively. 
Radial velocities (12) were measured from our spectra. 



3.2. A Ori and Orion 0B1b 

As in the case of ONC, we selected stars of A Ori from a wider 
sample acquired to investigate membership and accretion prop- 
erties in this young cluster by Sacco et al. (2008). To select 
the stars, we applied the same criteria used for the ONC sam- 
ple, discarding the fast rotating members (vsin/ > 30 km/s), 
binaries, and stars with spectral-type later than K7; this crite- 



rion resulted in the selection of one star only (Dolan 24). The 
observations were caiTied out in Service Mode during October 
and November, 2005, using FLAMES-UVES and CD#3 cross- 
disperser, covering the spectral interval ~ 4840 - 6840 A. We 
enlarged our sample with a member of the OB lb association 
(HD 294297) whose FLAMES-UVES spectrum, obtained with 
the CD#4 cross-disperser, was retrieved from the ESO Archive 
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Fig. 2. Same as Fig.[T] but UVES spectra of the A Ori star (upper panel) and of the OB lb member are shown. Note the different 
spectral interval displayed here with respect to Fig.[T] 



Table 2. A Ori and OB lb targets, observing log, and target properties' 



Star Obs. Dates Instrument Exp. time S/N I R-I Sp.Type Teff Vrad 

(sec) (K) (km/s) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



HD 294297 2005-11-17 



UVES, CD4 3x3575 150 



GO 6100 26.46 + 0.41 
— 5080 22.83 + 0.53 



Dolan24 2005-10-16/17/18 UVES, CD3 4x2760 200 11.75 0.53 — 

2005-11-11 

Q Information for star HD 294297 comes from the SIMBAD database. Identifiers and photometry of Dolan 24 in A Ori were taken from Dolan 
& Mathieu {T999l l. Teff for HD 294297 was retrieved from Cunha et al. ( fT998t , while that of Dolan 24 was estimated from the R-I color. Radial 
velocities were measured from our spectra. 



(Program 076.C-0145, PI Jeffries). This star is included in the 
sample of Cunha et al. ( 1998) and thus allows a direct compari- 
son of their and our abundance scales. For both stars data reduc- 
tion was performed under ESO-MIDAS context and using the 
FLAMES/UVES pipeline. The two targets along with informa- 
tion on their properties and observations are listed in Table ID 
their spectra are shown in Fig.|2l 

4. Analysis and abundance determination 

4.1. Radial velocities 

Radial velocities were computed from our spectra using the task 
RVIDLINES in IRAFQ and a list of 10 lines; with the exception 



of Dolan 24, all measurements are in good agreement with previ- 
ous determinations in the literature (Sicilia-Aguilar et al. 120051 
Caballero I2007I I. For Dolan 24 we instead obtain a radial ve- 
locity lower than the measurement of Dolan & Mathieu (119991 
-29.14 km/s) and slightly below two different estimates of the 
average radial velocity of A Ori (Vrad = 24.3 + 2.8 km/s -Dolan 
& Mathieu [T9991 v^ad = 27.3 + 0.49 km/s -Sacco et al. 2008|, 
suggesting that this star likely is a spectroscopic binary. Derived 
radial velocities for all stars are listed in the last columns of 
Tables[I]and|2] 



^ IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of Universities 



for Research in Astronomy, under contract with the National Science 
Foundation. 
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4.2. Estimate of veiling 

As well known, spectra of PMS stars might be affected by spec- 
tral veiling that alters measured equivalent widths (EWs), since 
photospheric absorption lines are filled in by accretion shock 
emission. As a consequence, the observed lines are weaker than 
intrinsic ones and one needs to correct for veiling before deriving 
abundances. This is accomplished by determining r, the ratio of 
the excess to the photospheric continuum; then the relationship 
between the true and measured EWs is EWo = EWmeas. x (1 H- r). 

As mentioned in the previous section, veiling for the two 
ONC stars observed with Giraffe was already derived by Palla 
et al. (| 2005 ll; both of them have r 0, and this was indeed the 
reason why we included them in our sample. On the other hand, 
we directly estimated r from our spectra for the seven ONC stars 
observed with UVES and for the A Ori member (HD 294297 
should not be affected by veiling -Cunha et al. 1995), following 
a procedure similar to that described in Palla et al. (2 005. 2007) . 
but using different features. Namely, we measured the EWs of 
selected strong lines in the target stars and compared them with 
those measured in the spectra of 1 1 members of the IC 2602 and 
IC 2391 clusters, which are old enough (~30-50 Myr) to ensure 
that their spectra are not affected by veiling. The IC clusters stars 
cover the same range of effective temperatures as the Orion tar- 
gets and their spectra are characterized by a similar resolution 
(see Stauffer et al. [19971 Randich et al. l2U0T\ . For the 580 set- 
up we used nine lines; namely, Ca i 5857.5 A, Ca i 6102.7 A, 
Ca I 6122.2 A, Fe i 6546.3 A, Ni i 6643.6 A, Fe i 6662.5 AA, 
V I 6624.8 A, Ca i 6717.7 A, Ti i 6743.2 A. For the 860 nm set- 
up we employed six lines in the range 6743-8085 A; namely, 
6743. 12 A (Fe i), 7 148. 15 A (Ca i), 7445.75 A (Fe i), 7698.97 A 
(K i), 7937. 13 A (Fe i), 7998.97 A (Fe i). 

As a first step, we checked whether a dependence of the 
strength of the lines on temperature was present among IC clus- 
ter members: with the exception of the Ca i 6122.2 A and K I 
7698.97 A features, all the other lines did not show strong trends 
with Teft within ~ 300 K. Hence, in these cases we derived the 
mean EW considering all stars in the two clusters and used this 
average value for comparison with the Orion stars. For the two 
features whose strength depends on effective temperature, we in- 
stead used the EW value of the IC cluster star with Teft closest to 
that of the given ONC star. For each line, the quantity nine was 
calculated as E Wic /E Wqnc - 1 ■ We did not find any trend of nine 
with wavelength for none of the stars, similarly to the results of 
Santos et al. (2008 ). Hence, we computed the average value of r 
for each star using all lines. These final values are listed in Col. 2 
of Tables [3] and |4] The two tables show that r is consistent with 
zero for all stars with exception of the ONC members h278, for 
which we obtained r - 0.8 + 0.07; due to this relatively high 
veiling, the star was discarded from the sample. 

4.3. Abundance analysis 

The analysis was performed using MOOG (Sneden fl 9731 -2002 
version) and a grid of 1-D model atmospheres from Kurucz 
( II993I I, by means of both EWs and spectral synthesis. Radiative 
and Stark broadening are treated in a standard way in MOOG; 
as for collisional broadening, we used the Unsold approxima- 
tion (119551 1 for all the Unes. As discussed by Paulson et al. 
( 1200 3') this choice should not greatly affect the differential anal- 
ysis with respect to the Sun. We also mention that very strong 
lines (EW> 150 mA) that are most affected by the treatment of 
damping have been excluded from our analysis. 



4.3.1 . EW analysis and stellar parameters 

As mentioned earlier, the spectrum of the A Ori star was ac- 
quired using the CD3 and the analysis was done using the line 
list of Randich et al. (2006) who had used the same UVES set- 
up. For the ONC targets and the spectrum retrieved from the 
archive that were observed with the CD4, we instead built a 
new line list. Although the spectra cover the wavelength inter- 
val 6700 to 10000 A, only the blue part up to ~ 8000 A was 
usable, since the red part was contaminated by several telluric 
lines which could not be corrected, owing to the lack of spec- 
tra of early-type stars. Fe i features to be used for the analysis 
in the interval ~ 6700 - 8000 A were selected from different 
sources in the literature and subsequently checked for suitability 
(e.g. for blends) both on the solar spectrum observed with UVES 
and on the best quality sample spectra. The final list contains 40 
Fe I lines for the warm star HD 294297 and 26 lines for the cool 
members of the ONC (see Table |5]l; 24 lines are in common. 
Not for all the cool stars, however, it was possible to measure all 
the 26 lines. Most log gf values were taken from Clementini 
et al. ( 119991 1; for the few lines not included in that study, 
log gf were instead retrieved from the Vienna Astronomical 
Line Database (VALE|3). We performed the analysis of the so- 
lar spectrum obtained with UVES at the same resolution as that 
of our sample spectra using both the line list with 40 lines and 
that with 26 lines. Adopting as solar parameters Tgff q = 5770 K, 
log go = 4.44 and = 1.1 km/s (see Randich et al. 2006), with 
both line lists we obtained log n(Fe)=7. 50+0.03, a value very 
close to log n(Fe)0=7.52 (Anders & Grevesse, 1989). Individual 
logn(Fe) values for the Sun are listed in Table |5] [Fe/H] for the 
sample stars with CD4 spectra was derived differentially with 
respect to our own determination. For Dolan 24, analyzed using 
the line list of Randich et al., a value logn (Fe)© = 7.52 was 
instead adopted. 

EWs were measured using the package SPECTRE and a 
gaussian fitting procedure. Although the spectra had been previ- 
ously normalized, local continuum was inspected and, if needed, 
adjusted at each EW measurement. Initial stellar parameters 
were estimated as follows: /. effective temperatures for the ONC 
stars were retrieved from Hillenbrand (1997 ) who, in turn, had 
derived them from spectral-types using the scale of Kenyon & 
Hartma nn ([1995l l. Teff of HD 294297 was taken from Cunha 
et al. (II9981 I, while Teff for the A Ori star, for which the 
spectral-type is not available, was derived from the R-I color 
and again the scale of Kenyon & Hartmann. A reddening E(R- 
I)=0.07 was adopted (see discussion in Dolan & Mathieu 1999^; 
//. An initial microturbulence ^=1.5 km/sec was adopted for 
all stars; Finally, surface gravities for the ONC stars and 
Dolan 24 were estimated from the relation between M, L, and 
Teff (logg=4.44+log(M/M0)-log(L/Lo)+41ogTeff-15.O447). 
For the ONC stars masses and luminosities were taken from 
Getman et al. ( 120051 ); for the star Dolan 24 in A Ori we derived 
the bolometric luminosity from the Ic magnitude, applying the 
bolometric correction by Bessell (1991) and adopting a distance 
to the cluster equal to 400 pc. Mass was then estimated using the 
PMS evolutionary tracks by Baraffe et al. (1 19981 1 with a = 1 . For 
HD 294297 we adopted the log g value of Cunha et al. ( fT998l l. 

As usually done, final temperatures were derived by remov- 
ing the trends between abundance values (log n(Fe)) and exci- 
tation potential (x) in MOOG. Similarly, final ^ values were de- 
rived by removing the trend between log n(Fe) and measured 
EWs. In both cases a 2cr clipping was applied to the initial line 



^ http://ams.astro.univie.ac.at/cgi-bin/vald/ 
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list before removing the trends. On the other hand, because of the 
lack of suitable Fe ii lines, log g values could not be optimized 
for any of the stars and initial values were adopted. 

Final stellar parameters are listed in Tables |3] and |4] We first 
note that for HD 294297 we retrieve the same Tefr as Cunha et 
al. ( I1998I I. Then, the comparison of Tables[T]and|2]with Tables[3] 
and|4]shows that the maximum difference between final and ini- 
tial effective temperatures is on the order of 100 K, with the ex- 
ception of star h673 (ATeff ~ 300 K) and Dolan24 (ATeff =270 K). 
All final temperatures are warmer than initial ones. Tables[3]and 
|4] also indicate that the final microturbulence values are in the 
range ~ 1.5-1.7 km/s, only slightly larger than those normally 
found for main sequence dwarfs of similar temperature. For the 
PMS stars in their sample, and in particular for the three ONC 
members, Santos et al. (?2008) instead found higher values (up 
to ^ = 2.5 km/sec) and a larger star-to-star scatter. They sug- 
gested that the discrepancy between their microturbulence val- 
ues and those of main sequence stars might be due to the effect 
of strong magnetic fields that characterize very young stars. Our 
results do not seem to support this hypothesis, since the micro- 
turbulence values for our sample stars are only slightly higher 
than that of the main sequence member of IC 2391 (see Sect. 14.41 
below). Note that use of higher values of the microturbulence 
would lead to lower [Fe/H] values. 

The iron abundance for each line was obtained based on 
measured EWs and adopted stellar parameters. Final abundances 
for each star were then determined as the mean abundance from 
the different lines. 

4.3.2. Spectral synthesis 

The spectral synthesis for the Fe analysis was performed in a 
20 A wavelength range around the Li i feature at 6707.78 A, 
from 6695 Ato 6715 A. As in the case of the EWs, we first car- 
ried out an analysis of the solar spectra acquired with UVES 
and Giraffe. When necessary, we modified \og gf values to get a 
good fit of the solar spectrum with the standard solar abundances 
of Anders & Grevesse ( I1989I ). 

The synthesis was carried out for all but one star observed 
with UVES, in order to cross-check [Fe/H] from EWs; spectral 
synthesis was also performed for the two ONC Giraffe spectra, 
whose lower resolution and reduced spectral range did not allow 
us to carry out the EW analysis. Also note that, due to the rela- 
tively high rotation, [Fe/H] for h907 in ONC could be measured 
only from the synthesis analysis. Instead, for the star h223a in 
the ONC we could not apply the synthesis method because part 
of the bluest order, including the wavelength range from 6695 to 
6740 A, fell out of the CCD. 

For stars observed with UVES, the synthesis was carried out 
adopting the final stellar parameters determined with the EW 
analysis; on the other hand, stellar parameters could not be opti- 
mized for stars observed with Giraffe and for h907, for which we 
did not carry out the EW analysis. Therefore, for these stars we 
adopted initial Tefr values and a microturbulence ^=1.6 km/s, 
which is typical of our UVES sample. Also, for each UVES 
spectrum, we first computed a synthetic spectrum with the metal- 
licity determined through the EW analysis; then, if necessary, we 
changed the metallicity until the best fit of the observed spectrum 
was obtained. For the analysis of the Giraffe spectra and that 
of h907 we started with a solar metallicity. Examples of spec- 
tral synthesis in the 6695-6715 A region for UVES and Giraffe 
spectra are shown in Figs.[3](star h673),|4](star h683), and|5](star 
h664), respectively. In all figures, both the whole spectral region 



and a zoom around two Fe i lines are displayed. Those two lines 
have different excitation potentials and thus their strength has a 
different dependence on effective temperature. A large error in 
the determination of the latter would thus be identified in the 
comparison of synthetic and observed spectra. We also note that 
the red side of the Li i 6707.8 A feature is not well fitted in the 
spectra of stars h683 and h664 (Figs. |4] and |5]). Although a de- 
tailed investigation of this effect is outside the purposes of this 
paper, we suggest that it might be due to the contribution of ^Li 
to the lithium feature. ''Li is indeed not included in our current 
line list. We finally notice that, as a by-product, spectral syn- 
thesis also allowed us to estimate projected rotational velocities 
(vsin /) for the sample stars. The latter are listed in Column 6 of 
TablesSandgl 

4.4. Errors 

Three sources of random errors are present in the EW analysis: 
(i) cTi, the standard deviation around the mean due to differences 
in abundances that come from each line; cri should be repre- 
sentative of errors in EWs and log gf values; (ii) 1x2: errors due 
to uncertainties on the stellar parameters. We found that for all 
the stars variations in Tgff greater than 100 K would produce 
significant trends in log n(Fe) versus excitation potential, while 
variations greater than 0.15 km/s in ^ would introduce a trend 
between log n(Fe) and measured EWs; we thus assumed 100 K 
and 0. 15 km/s as typical uncertainties in Teff and ^. As to gravity, 
we adopted a typical uncertainty of 0.25 dex. Resulting errors in 
[Fe/H] were then estimated by varying one parameter at the time 
and quadratically adding the resulting errors. We found that a 
change in Tgff of + 100 K results in a change of Fe abundance 
from +0.01 to +0.04 dex; a change in ^ of + 0.15 km/s causes a 
change of ±0.03 dex, and a variation of 0.25 dex in log g results 
in a change from ±0.04 to ±0.07 dex in abundance; (iii) 0-3: er- 
rors due to uncertainties on veiling determination. As shown in 
Tables [3] and H) the last are very small and lead to an error in 
the final [Fe/H] on the order of 0.01 - 0.02 dex, much below the 
other two errors. Also note that those small errors on veiling do 
not introduce any effect on the determination of stellar parame- 
ters. 

The abundance uncertainties are related to the best fit de- 
termination (crl) and to the uncertainties in stellar parameters 
(cr2), also in the case of spectral synthesis analysis. Assuming 
the same uncertainties in stellar parameters derived from the EW 
analysis, we found that a change in Teff of ±100 K results in 
a difference in Fe abundance from ±0.05 to ±0.1 dex, while a 
change of ^ of ± 0.15 km/s causes a variation from ±0.03 to 
±0.05 dex in abundance values. The effect of gravity is instead 
negligible (below 0.02 dex). 

Systematic errors that could affect the analysis of cool stars 
were estimated by analyzing with the same method a star with 
similar parameters as our sample stars and known metallicity. 
We chose a member of the open cluster IC 2391 (VXR76A) 
whose metallicity has been estimated with an independent 
method (different line list and code) by Randich et al. ( 20011 - 
[Fe/H]= 0.0 ± 0.09 ± 0.09). For the analysis of VXR76A we 
used the same spectrum employed by Randich et al., charac- 
terized by a spectral resolution R=43800, very similar to that 
of the UVES spectra. We derived the metallicity of VXR76A 
using both EWs and the two line lists and spectral synthesis. 
The results are listed in the last raw of Table [3] In all cases the 
agreement with the previous determination is excellent: namely, 
[Fe/H]= 0.0 ± 0.06 ± 0.08 from the EWs and the CD3 line list. 
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Fig. 4. Same as Fig. [3] but the observed and synthetic spectra of star h683 are shown. Note that the [Fe/H] from EWs for this star is 
slightly larger than that from the synthesis; the figure shows that the synthetic spectrum with [Fe/H] =-0.1 better fits the observed 
spectrum. 



[Fe/H]= -0.01 +0.08+0.08 from the EWs and CD4 Hne Ust, and 
[Fe/H]= + 0.1 +0.1 from the spectral synthesis. Furthermore, 
in a more recent analysis of a wider sample of stars and better 
quality spectra of IC 2391, D'Orazi & Randich (2009) not only 
confirmed a solar metallicity for this cluster, but also found that 
VXR76A shares the same metallicity as other, warmer members 
in IC 239 1 . This suggests that our analysis should not be affected 
by major systematic eiTors and that no major offset should be 
present between the metallicity of cool and warmer stars. As for 
Dolan 24, the method/line list was already tested by Randich et 
al. (I2006I I on M 67 members. 



5. Results and discussion 

5.1. The metallicity of the ONC 

The results of our analysis are summarized in Tables |3] and |4] 
where we give the veiling value r (Col 2), the adopted Teff , log g. 



and f values (Cols. 3-5), vsin/ (Col. 6), [Fe/H] derived from 
EWs analysis (Col. 7), along with the number of lines employed 
for the analysis. [Fe/H] values from spectral synthesis are listed 
in Col. 8. 

First, the two tables show that, in general, [Fe/H] values 
from EWs and spectral synthesis well agree within the uncer- 
tainties. We derived the mean ONC metallicity considering the 
more secure [Fe/H] from EWs for stars with both determina- 
tions; [Fe/H] from synthesis was instead considered in the three 
cases for which the EW analysis could not be performed. We 
find [Fe/H]av = -0.01 + 0.04, i.e., the ONC has a solar metalUc- 
ity. Including only stars with [Fe/H] from EWs, we would obtain 
[Fe/H]av = -0.02 + 0.05. 

As mentioned, Santos et al. (I2008I I. based on three stars, 
derived a significantly lower average metallicity for the ONC 
([Fe/H]= -0.13 + 0.06), although the [Fe/H] of the most metal- 
rich star in their sample is comparable to the most metal- 
poor one in our own sample. The discrepancy is likely due to 
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Fig. 5. Same as Figs. |3] and |4] but the observed and synthetic spectra of star h664 observed with Giraffe are shown. Note that Giraffe 
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Table 3. Final stellar parameters and metallicity values derived for ONC stars. ' 



star 


r 


Teff 


logg 




vsin i 




[Fe/H] + 0-1 


+ (T2 






(K) 




(km/s) 


(km/s) 




(EWs) 


(Spec. Synth.) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


223a 


0.02 ± 0.04 


4450 


3.8 


1.6 




0.03 ± 0.06 


+ 0.07 (18) 




268 


0.03 ± 0.04 


4300 


3.9 


1.6 


10+3 


-0.08 + 0.16 


+ 0.09 (25) 


-0.1+0.1+0.1 


278 


0.8 ± 0.07 


4395 














487 


0.02 ± 0.02 


4300 


3.9 


1.7 


15±3 


-0.07 ± 0.08 


+ 0.09 (18) 


+ 0.1+0.09 


673 


0.04 ± 0.03 


4700 


4.0 


1.5* 


20+3 


0.0 + 0.1 


+ 0.07 (18) 


+ 0.1+0.08 


683 


0.0 ± 0.02 


4250 


3.3 


1.6 


13±3 


0.01+0.1 


+ 0.08 (22) 


-0.1+0.1+0.09 


907 


0.01 ± 0.03 


4775 


4.1 


1.5 


35±3 






0.0 + 0.1+0.11 


664 


0" 


4200 


4.1 


1.6 


13±3 






0.05 + 0.1+0.1 


1020 




4200 


3.8 


1.6 


13±3 






0.0 + 0.1+0.1 


VXR76A (listcDs) 




4340 


4.5 


1.2 


8±3 


0.00 + 0.06 


+ 0.09 (25) 


+ 0.1+0.1 


VXR76A (listcD4) 




4400 


4.5 


1.3 


8±3 


-0.01+0.08 


+ 0.08 (23) 


+ 0.1+0.1 



'The initial microturbulence was adopted for this star, since £ could not be optimized due to the lack of 

lines covering a wide enough EW interval. 

"r values for stars 664 and 1020 come from Palla et al. (2005 1. 
Un the table we list the star ID, the veiling r, adopted stellar parameters, projected rotational velocities, and the final [Fe/H] values, from the EW 
and synthesis analysis, along with cri and o"2 errors (see text). For [Fe/H] from EWs we also provide in parenthesis the number of employed Fe I 
lines. 

Table 4. Same as Table|3] but the results for A Ori targets and star HD 294297 are listed. 



star 


r 


Teff 


logg 




vsin I 






[Fe/H] + o-i 


+ 0-2 






(K) 




(km/s) 


(km/s) 






(EWs) 


(Spec. Synth.) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 






(V) 


(8) 


HD 294297 




6100 


4.0 


14 


<7.5 


-0.16 


+ 0.03 


+ 0.07 (34) 


-0.15 + 0.05 + 0.08 


Dolan24 


0.05+0.03 


5350 


4.3 


1.4 


<7.5 


0.01 


+ 0.02 


+ 0.06 (54) 


0.1+0.1+0.08 



the higher microturbulence values obtained in that study (see 
Sect. 14.3.1b . As in our case, Santos et al. have used FLAMES 
with the fiber link to UVES and their spectra are of similar 
quality. For the analysis they have employed the same code and 
model atmospheres; however, they adopted a different line list 
and cr-clipping criterion. We suggest that this affects the mi- 
croturbulence determination and thus the final [Fe/H]. We stress 
again that the comparison with VXR76A provides confidence in 



our results. As discussed in Sect. 1, the metallicity of the ONC 
has so far been poorly constrained, with some studies suggest- 
ing an underabundance with respect to the Sun and others finding 
solar or even over-solar [Fe/H] values. Iron has been determined 
in eight low-mass stars (1 from Cunha et al. 119981 four from 
Padgett [T996l and three from Santos et al. 120081 ) and one B-type 
star (Cunha & Lambert 1994 ), but the [Fe/H] values are not ho- 
mogeneous, since they were measured with different techniques. 
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Our study, not only almost doubles the number of ONC low- 
mass members with a metallicity measurement, but our [Fe/H] 
are homogeneous, thus allowing us to secure the average metal- 
licity of the ONC on more solid grounds. The left panel of Fig.|6] 
shows the distribution of [Fe/H] for the sample of Cunha et al. 
( 119941 119981 ) including both early- and late-type stars in all the 
four subgroups; in the middle panel we plot the [Fe/H] distribu- 
tion of ONC only, considering all previous determinations from 
the literature; finally, the distribution of our own measurements 
for ONC members is displayed in the right-hand panel, along 
with our average and the typical error on [Fe/H] of individual 
stars. The figure shows that the distribution is very broad when 
including other subgroups besides the ONC, possibly implying 
a group-to-group difference; the distribution becomes narrower 
when considering the ONC only and measurements from differ- 
ent groups, but still spans ~ 0.4 dex in [Fe/H]. The distribution 
of our own [Fe/H] values is instead much narrower and the small 
dispersion still present is completely due to measurement uncer- 
tainties. 

As to OB lb, our [Fe/H] determination for star HD 2 94297 
is in excellent agreement with the value of Cunha et al. ( 119981 
-[Fe/H]= -0.19). For the same star Gonzalez-Hernandez et 
al. ( I20U81 ) derived [Fe/H]= -0.09 + 0.1 (star SO000041 in their 
list), also in good agreement. All the three measurements in- 
dicate a sub-solar metallicity, much below our average for the 
ONC. In our study the [Fe/H] values for this star and the ONC 
are on the same scale; we mention in passing that, by using for 
the analysis of HD 294297 the same subset of lines used for 
the ONC, we would still obtain for this star the same [Fe/H]. 
This suggests that the difference between the ONC stars and 
HD 294297 is real and that the OB lb group might be more metal 
poor than the ONC. We further develop this point in the next sec- 
tion. 

Finally, the metallicity of A On is likely solar; whereas this 
represents the first measurement of the metallicity in this cluster, 
it is based on one star only and, obviously, additional measure- 
ments are needed. 

5.2. The iron distribution in Orion 

In Fig. |7] we plot the spatial distribution of stars in the differ- 
ent Orion subgroups with an available metallicity determination; 
stars in three metallicity bins are indicated with symbols of dif- 
ferent size. For the ONC we show our average measurement. 

The figure indicates that the three stars belonging to the 
OB lb subgroup all have [Fe/H] below the average of the ONC 
and, possibly, of the Ic subgroup, whose metallicity is charac- 
terized by a larger uncertainty and dispersion. Two of the stars 
in the OB lb sample are B-type stars from Cunha & Lambert 
d 1 9941 1, while the third one is HD 294297; for this star we adopt 
our own [Fe/H] determination. The resulting average metallic- 
ity for the lb subgroup is [Fe/H]= -0.10 + 0.05, 2cr below 
our determination for the ONC. The difference is not big and 
based on very few objects only, two of which are early-type 
stars. Nevertheless, it suggests that the Orion region may not 
be characterized by a single value of [Fe/H]. As well known, 
Fe is mostly produced by type-I Supernovae (SNe), whose life- 
times are much greater than the age of the Orion association 
(~ 10^ years). On the other hand, massive, short-lived stars ex- 
ploding as type-II SNe might pollute a molecular cloud close to 
OB associations. Although SNe II are mainly producer of a el- 
ements like O, Mg, and Si, a small amount of Fe is also pro- 
duced. As mentioned in Sect. |2l Cunha & Lambert ( 119921 ) and 
Cunha et al. dl 9981 1 reported evidence of star-to-star variations 



in O and Si abundances within Orion, with a few O/Si-enhanced 
stars more centrally located in the Trapezium region than the 
widespread low O/Si ones; they suggested a simple scenario of 
supernova self-enrichment across the Orion association over the 
past 10 Myr. With the caveat that the Fe yields from SNe II are 
highly uncertain, they estimated that the enhancement in oxy- 
gen of ~ 0.3 dex would correspond to an enhancement in Fe of 
~ 0.06 dex. Whereas in their study they did not find such a vari- 
ation, the ~ 0.1 dex difference in iron that we measure between 
the ONC and lb subgroup is consistent with this the scenario. 
The OB 1 a subgroup might also be metal-poor; three of the four 
members (all from Cunha & Lambert |1994| l have [Fe/H] more 
than Icr below the ONC, while one star is overabundant with 
respect to the Sun and the ONC. 

Finally, Fig. Q shows that a dispersion in metallicity is 
present in OBlc. We suggest that the dispersion is likely due 
to both systematic differences between different studies (i.e., 
Cunha et al. versus Padgett) and to large uncertainties within 
a given study, as we have shown to be the case for the ONC. 

5.3. Tlie lack of metal-ricli Stars 

One of the goals of our project is the search for metal-rich young 
low-mass stars. The first results presented here indicate that the 
ONC has a solar metallicity; the dispersion we find within the 
cluster is consistent with measurement errors and none of the 
stars in our sample is metal-rich. Under the assumption that our 
analysis is not affected by major systematic errors, as the com- 
parison with IC 2391 confirms, two opposite hypothesis can be 
suggested to explain the lack of such a population. The simplest 
and most straightforward one is that they do not exist and, as ex- 
pected within a cluster, the ONC is characterized by a homoge- 
neous composition of all its members. Alternatively, the [Fe/H] 
distribution of the ONC is inhomogeneous and metal-rich mem- 
bers exist, but low number statistics prevented us from detecting 
them. To this respect, we recall that one star from Padgett d 1 9961 ) 
has a metallicity significantly above solar ([Fe/H]=0.I4 + 0.18) 
and above our mean for the ONC. The final answer on the possi- 
ble presence of metal-rich stars can only be provided by a metal- 
licity determination in a much wider sample of stars, including 
also a re-analysis of the possible metal-rich ONC member from 
the study of Padgett. We caveat however that a large metallicity 
dispersion within the ONC itself would be difficult to explain 
without assuming ad-hoc processes. To our knowledge, abun- 
dance patterns in all open clusters so far investigated are rather 
homogeneous and no metal-rich members of solar-metallicity 
clusters or, viceversa, solar-metallicity members of metal-rich 
clusters have been detected. 

On more general grounds, the question arises on the ex- 
istence of metal-rich young associations and SFRs. None has 
been detected so-far (see the discussion in Santos et al. 1 20081 ). 
In Fig. |8]we show the metallicity distribution of open clusters 
within 500 pc from the Sun with an available spectroscopic 
[Fe/H] determination. The 16 clusters included in the sample 
have ages in the interval 30 Myr-2 Gyr; the metallicity for all 
but one has been derived from high-resolution spectra of low- 
mass stars and typical errors on the average [Fe/H] values are 
in the range 0.02-0. 1 dex. The figure shows that the distribution 
peaks at solar [Fe/H], with a not negligible tail at high values; 
specifically, metallicities vary between [Fe/H]~ -0.15 and 0.4 
with an average value [Fe/H] = 0.02 + 0.1. In other words, the 
majority of the clusters in the solar neighborhood share the same 
metallicity as the Sun and the ONC is not peculiar. On the other 
hand, metal-rich open clusters do exist in the solar circle, and 




Fig. 6. Left panel: [Fe/H] distribution across all the four Orion subgroups from Cunha et al. (19941 [19981 ); middle panel: [Fe/H] 
distribution of the ONC only, obtained considering four different studies (Cunha & Lambert 1994; Cunha et al. J_998; Padgett 1996J 
Santos et al. 2008 ); right panel: [Fe/H] distribution of the ONC from our study. We also indicate our average for the ONC (vertical 
line) and the typical random error on [Fe/H] of individual stars. 



their number, and thus probability of finding them, is in principle 
not negUgible (3/16 or 19 % of the clusters have [Fe/H]> 0.1). 
However, whereas a strict [Fe/H] versus age relationship does 
not hold, the four clusters with [Fe/H]> 0.05 are all older than 
200 Myr; furthermore, two of them (Hyades and Praesepe) be- 
long to the Hyades super-cluster or stream (Eggen [T992[ Famaey 
et al. I2007I I and might have a peculiar origin. The majority 
of the clusters with close-to-solar metallicity are instead part 
of the Local Association (or Pleiades moving group) or of the 
IC 2391 supercluster, to which various young associations such 
as Scorpius-Centaurus are also associated. Finally, the open clus- 
ters IC 4665 and IC 2391, both with a solar iron content are 
likely members of the Gould Belt (Piskunov et al. 120061. In 
other words, both kinematic and age considerations suggest that 
it might not be easy to find a metal-rich young association in the 
solar neighborhood. 



larger than 500 pc might indeed be feasible with current instru- 
mentation. 
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6. Concluding remarks 

Using FLAMES/UVES and Giraffe spectra, we determined iron 
abundances in eight low-mass members of the ONC, one star in 
the A Ori cluster, and one member of the Orion OB lb subgroup. 
The average metallicity of the ONC is [Fe/H]= -0.01+0.04; the 
metallicity of Dolan 24 in A Ori is also solar, while the star be- 
longing to the OB lb subgroup is significantly more metal-poor 
than the ONC. The ~ 0.1 dex difference in the Fe content of 
the ONC and the OB lb subgroup is consistent with the self- 
pollution scenario originally proposed by Cunha and collabo- 
rators ( fT992l [19941 IT998] ) for the Orion region. 

None of our sample stars nor any of the SFRs so far studied 
is metal-rich. At the same time, known metal-rich open clusters 
are all older than 200 Myr, and none is spatially or kinematically 
associated to young nearby associations and SFRs. This suggests 
that finding the metal-rich SFRs among the best studied ones in 
the solar neighborhood might not be easy or doomed to failure. 
We suggest that in order to find them, one should possibly look 
outside the Gould Belt and the Local Association. High-spectral 
resolution observations of low-mass stars in SFRs at distances 
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Fig. 7. Spatial distribution of members of the Orion association with [Fe/H] measurements. Open symbols indicate members of the 
OBla,lb, Ic subgroups, while our average for the ONC is represented by the star symbol. Symbols with different size denote three 
different metallicity bins: [Fe/H]< -0.05, -0.05 <[Fe/H]< 0.03, and [Fe/H]> 0.03, from the smallest to the largest bin. The centi-al 
bin corresponds to the average metallicity of the ONC +lcr. The dashed line outhnes the boundaries of the Ori OB lb subgroup 
following Warren & Hesser (|1977l l. 
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Fig. 8. Metallicity distribution of clusters within 500 pc from the Sun with a secure measurement of the metallicity; the average ±lo- 
for the ONC is indicated as a vertical bar. 
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Table 5. Line list of 42 Fe i lines used for the analysis of the ONC stars along with excitation potential log gf values, equivalent widths 
measured on the solar spectrum acquired with UVES and corresponding Fe abundances. Lines labeled with a 'w' or 'c' have been used for the 
warm star HD 294297 or for the cool ONC stars only. 
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